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ABSTRACT 

We have developed a new method to extract halo merger rates from the Millennium Simulation. First, by 
removing superfluous mergers that are artifacts of the standard Friends-Of-Friends (EOF) halo identification 
algorithm, we find a lower merger rate compared to previous work. The reductions are more significant at 
lower redshifts and lower halo masses, and especially for minor mergers. Our new approach results in a better 
agreement with predictions from the extended Press-Schechter model. Second, we find that the EOF halo 
finder overestimates the halo mass by up to 50% for halos that are about to merge, which leads to an additional 
« 20% overestimate of the merger rate. Therefore, we define halo masses by including only particles that are 
gravitationally bound to their EOF groups. We provide new best-fitting parameters for a global formula to 
account for these improvements. In addition, we extract the merger rate per progenitor halo, as well as per 
descendant halo. The merger rate per progenitor halo is the quantity that should be related to observed galaxy 
merger fractions when they are measured via pair counting. At low mass/redshift the merger rate increases 
moderately with mass and steeply with redshift. At high enough mass/redshift (for the rarest halos with masses 
a few times the "knee" of the mass function) these trends break down, and the merger rate per progenitor halo 
decreases with mass and increases only moderately with redshift. Defining the merger rate per progenitor halo 
also allows us to quantify the rate at which halos are being accreted onto larger halos, in addition to the minor 
and major merger rates. We provide an analytic formula that converts any given merger rate per descendant 
halo into a merger rate per progenitor halo. Finally, we perform a direct comparison between observed merger 
fractions and the fraction of halos in the Millennium Simulation that have undergone a major merger during 
the recent dynamical friction time, and find a fair agreement, within the large uncertainties of the observations. 
Our new halo merger trees are available at http://www.mpe.mpg.de/ir/MillenniumMergerTrees/ 
Subject headings: cosmology: theory — dark matter — large-scale structure of the universe — galaxies: evo- 
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1. INTRODUCTION 



(e.g. jHernauist 198 ^, and transform gala ctic morphology 
(e.g.l Naab & BurkerfcoOStlBournaud et ani2005h . Many pa- 



in the current paradigm for structure formation, the cold 
dark matter mode l (White & ReesI 119781: iBlumenthal et alj 
11984 iDavis et alJ[T985i: iSpringel et al.l l2006'). galaxy merg- 
ers play an important role in galaxy formation and evolu- 
tion. ^Galax^;_jnergers_drive towards central starbursts 
(e.g.lMihos & Hern auistl[T996h and supermassive black holes 

ly 

pa- 
rameters, like the gas fraction and morphology of the merg- 
ing galaxies, as well as their relative orbits and orientations, 
affect the properties of the merger remnant. One of the most 
impor tant factors is t he mass ratio of the merging galaxies 
(e.g. Naab et al.ll2006l) . Galaxy merger s of all ma ss ratios are 
frequent in a ACDM Universe (e.g. Ste wart et an f2008b). but 
special importance is given to major mergers, usually consid- 
ered as those with mass ratios less than « 3 : 1. These are 
thought to play a significant role in the buil dup of the red se- 
quence jToomrdl 19771: iHopkins et al.ll2008l) by transforming 
blue star-forming late type galaxies into red passive early type 
galaxies. 
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Many observational studies have been carried out in re- 
cent years to investigate the fraction of galaxies that show 
signs of major merger activity as a function of mass, lu- 
minosity and redshift. Observationally, only merger frac- 
tions can be obtained, and in order to transform them into 
merger rates, which can be directly compared with theoretical 
models, th e time scale of the observed events must be esti- 
mated (e.g. Patton et alj|2000l) . Two principal approaches are 
used to o bserve galaxy merger fractions. One is pair count- 
ing (e. g. IPatton et all 119971: iUe Fevre et alj 120001: iBeU et alj 
I2006al: iRvan et alj 120081: iLin et alj 120081) . i.e. identifying 
galaxies separated from one another by less that typically w 
20kpc. This method probes the pre-merger stage a nd is there- 
fore a "progenitor galaxy" merger fraction ( De Propris et alj 
l2007h . The second approach is identification of mergers 
through morph ological signatures such as asymmetry and 
tidal tails (e.g. [Le Fevre et al.l 120001: IConselice etal] 120031: 
iLotz et al] |2008y This method aims at identifying merg- 
ers in their relatively late stages, and is therefore a "descen- 
dant galaxy" merger fraction. For each of the principal ap- 
proaches different authors use different methods, as well as 
different selection criteria, which create systematic effects 
that are not easily comparable. Moreover, the intrinsic un- 
cert ainties of eac h method, and the effect of cosmic variance 
(e.g. lConselice et al. 2008) contribute further to the large scat- 
te r between the obtained r esults. 

iDe Propris et alJ (l2007h and lConselic3 (l2006l) have noticed 
that the merger fractions obtained with the two principal 
methods should be carefully defined, and cannot be directly 
compared to each other, as they are actually different quanti- 
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ties. Both authors proposed (different) simple conversions be- 
tween the two quantities. Nevertheless, both disregarded the 
fact that the mass difference between the progenitors and the 
descendant is up to a factor of 2 (in equal-mass binary merg- 
ers). That may result in very large differences in the number 
densities of progenitors and of descendants. Since the mea- 
sured merger /racf/on implicitly includes the information of 
the number of objects, different number densities are implic- 
itly included in the merger rate or merger f raction per progen- 
itor galaxy ver sus per descendant galaxy. iBell et all (l2006bh 
and iLot z et al.l ("2008) realised this, and have taken the dif- 
ferent number densities into account when comparing merger 
fractions obtained with different methods. We show in this 
paper that the different definitions of the merger rate, per de- 
scendant or per progenitor halo, lead to significantly different 
results, particularly for major mergers, high redshift mergers, 
and high mass mergers. 

Galaxy mergers follow their dark matter halo mergers, but 
the connection between the two is not straightforward. When 
two dark matter halos merge, the orbital angular momentum 
is transferred into internal degrees of freedom, while the more 
concentrated galaxies are at first not much affected. The 
galaxies lose relative angular momentum due to dynamical 
friction, and t hey start merging as well when they are more 
tightly bound ( Barnes & H ernquist 1992). All of the baryonic 
physics involved in galaxy mergers makes it difficult theoreti- 
cally to quantify the galaxy merger rate reliably (P. F. Hopkins 
et al. 2009, in preparation). Nevertheless, a first step towards 
quantifying the galaxy merger rate would be to understand 
the dark matter halo merger rate, which can be studied more 
r obustly. 

iLacev & Cold d 19931) have estimated the halo merger 
rate based on the extended Press-Sc hechter (EPS) mode l 
(iPress&Schechted 119741 : iBond etalJll 991: Bower| 1199 Ih . 
iNeistein & Deke l' ('2008b') and 'Zhang et al.i (ilOOSaHiave re- 
cently constructed EPS-based approximations that are self- 
consistent, and Zhang et alJ (12008b) have done so for the el- 
lipsoidal collapse model. All these models differ by factors 
of a few tenths up to a few. N-body simulations have only 
recently been used to study the halo merger rate with a large 
dynamical rang e (D'Onghia et al. 2008; Stewart et al. 2008a; 
iFakhouri&Md 12008. hereafter FM08). FM08 have investi- 
gated this problem ba sed on the large Millennium Simulation 
(ISpringel et al.ll2005h . They have found that the dark matter 
halo merger rate has an almost universal form that can be sep- 
arated into its dependencies on mass ratio, descendant mass 
and redshift. Nevertheless, the analysis of N-body simulations 
is also subject to uncertainties, in particular in identifying ha- 
los. FM08 present 3 ways to analyse the simulation, which 
differ from one another by 25%, and by a mor e signifi- 
cant amount compared to the iLacey & Cold (Il993h analyti- 
cal approximation. All of these investigations have quantified 
the halo merger rate per descendant halo. Recent work has 
also used N-body simulations to study mergers of subhalos 
(lAngulo et al.ll2008l:IWetzel et al.ll2008l) . 

In this paper we extract the dark matter halo merger rate 
from the Millennium Simulation using a new method for iden- 
tifying halos and mergers, as well as a new definition of the 
merger rate. In ^2.11 we review the Millennium Simulation 
and its post-processing. In ^2.21 we describe how we create 
merger trees th at ar e free of artificial effects not considered 
previously. In ^2.3l we define the merger rate per descendant 
halo, as well as the merger rate per progenitor halo, which 
is related to observed galaxy merger fractions when they are 



measured via pair counting. In ^ we describe our results 
and in ^we compare them to previous work. In ^we cal- 
culate the halo merger fraction in the Millennium Simulation 
and compare to observations. In ^we discuss our results and 
their relevance to observations, and summarise. 

2. EXTRACTING THE MERGER RATE FROM THE MILLENNIUM 
SIMULATION 

2.1. The Millennium Simulation and its merger trees 

The Millennium Simulation is a cosmological N-body sim- 
ulation following 2160-' dark matter particles, each of mass 
8.6 X l(fhr^ Mq, in a box of 500/i"' Mpc on a side, with 64 
generated output times ("snapshots") from z= 127 to z = 0. 
The cosmology is set to ACDM with U.,„ = 0.25, ^Ia = 0.75, 
Qh = 0.045, h = 0.73, n = 1 and erg = 0.9, which we will adopt 
throughout this paper 

In the Millennium Simulation there are two levels of struc- 
ture id entification. First, the Friends-Of-Friends (EOF) algo- 
rithm ( Davis et al.llI985l; with a linking length b = 0.2 of the 
mean particle separation) creates a catalogue of EOF groups 
in every snapshot. In the limit of a large number of parti- 
cles, EOF groups enclose the particles within iso density con- 
tours of » b~^ times the mean matter density dErenk et alj 
IT988I: ILacev & Cold IT991 I.Tenkins et al.1 llOOlh . and their 
mean densities correspond approximately to the overdensi- 
ties of virialised halos expected from the spherical collapse 
model (Lacey & Cole 1994). Thus, EOF groups are consid- 
ered to represent dark matter halos. However, substructure as 
traced by local density maxima are not distinguishable within 
the EOF groups. For this purpose, the algorithm SUBEIND 
(ISpringel et al.l 1200 Ih identifies substructures in each EOF 
group, by finding gravitationally self-bound groups of parti- 
cles around maxima in their smoothed density field. Thus, 
each FOE group contains at least one subhalo, and the sub- 
halo wit h the "most massive history " is chosen to be the main 
subhalo dPe Lucia & Blaizot"2007^. 

The Millennium merger trees^ are constructed from the sub- 
halos by finding a single descendant for each subhalo at the 
following snapshot. The EOF groups themselves play no role 
in constructing the merger trees. Nevertheless, if the halo 
merger rate is to be studied, then new merger trees, in which 
each node is a halo rather than a subhalo, must be built. 

2.2. Constructing new halo merger trees 

We build new halo-based, that is EOF group-based, merger 
trees by defining one descendant for each EOF group. The 
main subhalo in each EOF group is identified and followed to 
its subhalo descendant using the original subhalo trees. Then, 
the EOF group to which that subhalo descendant belongs is 
defined as the EOF group descendant of the EOF group in 
question. 

In practice, EOF groups not only merge, but may also split. 
This is not compatible with the simplified notion of hierarchi- 
cal build-up and indeed not described at all by an analytical 
model such as EPS. Yet, this phenomenon is robust in numer- 
ical N-body simulations. A EOF group is described as "split" 
when a subhalo belonging to it at a given snapshot is no longer 
part of it at a later snapshot. 

Those subhalos that were split out of their EOF groups can 
be classified according to several criteria, describing their past 

' Structure catalogues and the derived merger trees have been made public 
by the Virgo Consortium: http://www.mpa-garching.mpg.de/millennium 
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Fig. 1. — Left: The distribution of relative distances (witli respect to S200 at each snapshot) between subhalos that only temporarily belong to the same FOF 
group. The sohd blue curve shows the distribution at the last snapshot of the temporary merger, the green curve with asterisks shows the minimum relative 
distances during the temporary merger, and in red with circles is the distribution of maximum distances. These distributions are shifted slightly towards lower 
(higher) values at lower (higher) redshifts and for more minor (major) mergers. Right: The distribution of the ratio between the fragmentation duration (the time 
between the fragmentation itself and the subsequent re-merger) and the first, temporary, merger that precedes the fragmentation, for cases where the fragmentation 
lasts more than three snapshots. The fragmentation duration is typically much longer (mean Ri 4, median ss 3) than the temporary merger duration, indicating 
that the temporary merger is usually only an artifact of the FOF algorithm. 



TABLE 1 

Merger scenarios, their abundance, and the different counting by different algorithms 



first merger duration in snapshots 


00 


< 3 


< 3 


< 3 


> 3 


> 3 


> 3 


fragmentation duration in snapshots 




< 3 


> 3 


00 


< 3 


> 3 


00 


abundance 


78% 


7.6% 


4.4% 


3.5% 


2.5% 


1.9% 


2.1% 


snipping 




2 


2 


1 


2 


2 




stitching-3 




1 


2 


1 


1 


2 




stitching- 00 




1 


1 


1 


1 


1 




splitting-3 




1 


1 





2 


2 




splitting 




1 


1 





1 


1 






Note. — In bold font ai'e scenarios for which a certain algorithm over-counts the number of mergers. For details see ^2.21 A merger duration of 00 snapshots 
denotes a merger with no fragmentation, and a fragmentation duration of 00 snapshots means that the fragments never re-merge by z = 0. The abundances of the 
different scenarios are calculated using all the mergers in the Millennium Simulation. 



and future histories: where they were created, what struc- 
ture they belong to right after the split, and where they end 
up at later times. Almost all subhalos are noi formed inside 
their host FOF groups, but originate from previous mergers in 
which they were accreted onto them. Thus, most split subha- 
los are part of mergers that started in the past but were cut by 
the split, before the subhalos merged. The location of the split 
subhalo in the immediate snapshot after the split can be either 
as an independent FOF group or as a subhalo in another FOF 
group. After the split, the split subhalo may re-merge with the 
FOF group from which it was split, or it may not do so at all. 

A very typical case is that of two FOF groups that merge, 
with one becoming a subhalo of the other, which then later 
split again into two distinct FOF groups. If they never re- 
merge, it is clear that a "merger" interpretation of such an 
event would not be appropriate. Most likely, a temporary 
bridge of particles caused the FO F algorithm to iden tify them 
as one FOF group (White 2001 iLukic et al.ll2008l) . If they 
eventually do re-merge after a few snapshots, the first merging 
and splitting merely shows that they were relatively close and 
that some interaction took place. However, the actual point 
of merging should correspond to the merger event after which 
there is no more splitting. To support that interpretation, we 
show in Figure [T(a)| that in 65% (85%) of temporary mergers 
the distance between the subhalos never decreases below .^200 
(7?2oo /2) of the joint FOF group. The undesired implications 



of the way splits appear in the original merger trees are artifi- 
cial changes in the halo mass function and merger counts. As 
shown by FM08 and as will be shown later on in this paper, 
different algorithms that deal with this phenomenon change 
the halo mass function only slightly, but have a pronounced 
effect on the merger rate. 

It is common not to treat such fragmentation cases in any 
special way, thereby allowing the split fragments to have no 
progenitors (e.g. FMOS's "snipping" method). In that case, 
the split fragment may re-merge (even several times) with 
the same halo. One approach of deali ng with this issue 
(e.g.' Hellv et"ai]|2003l:lHarker et al.ll2006l) is via different al- 
gorithms for splitting FOF groups that are artificial combi- 
nations of several halos. Another approach was taken by 
FM08, who chose in their fiducial "stitching-3" method to 
merge back the split fragments when a future merger takes 
place within three subsequent snapshots after the fragmen- 
tation. They also presented a method, "stitching-cx)", where 
subhalos are never allowed to split out of their FOF groups. 
In the present paper we eliminate the occurrence of halo frag- 
menta tion in our new trees, in a way similar to iHelly et aH 
( l2003h . We do so by identifying any FOF group at redshift 
Zp that contains at least a pair of subhalos that at some lower 
redshift z/ < Zp do not belong to the same FOF group. Such a 
FOF group is split by our algorithm into several fragments in 
the following way. Subhalos which belong to different FOF 
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groups at z/ will belong to different fragments at Zp as well, 
while subhalos that do stay together at z/ will not be sepa- 
rated at Zp- Any new fragment our algorithm creates, as well 
as any untouched FOF group, is considered hereafter simply 
as a "halo". IFakhouri & Mai (2009) (hereafter FM09) have 
suggested a variant of this algorithm, "splitting-3", that splits 
the progenitors of fragmenting FOF groups only 3 snapshots 
backwards. 

We identify two possible advantages of our "splitting" 
method over "snipping", "stitching" and "splitting-3" algo- 
rithms. The first has to do with double counting of merg- 
ers. Table [T] presents the possible scenarios for binary merg- 
ing and fragmentation events, and the number of mergers that 
each of the methods presented above counts for those scenar- 
ios. Printed in a bold font are cases where a method counts 
more mergers than appropriate. Only our "splitting" method 
never counts artificial mergers as real ones. Moreover, all 
other methods except from "stitching-oo" leave spurious frag- 
mentations in the tree, which may be mistakenly interpreted 
as extremely high (positive and negative) "smooth" changes in 
the halo mass, not associated with any merger. We note, how- 
ever, that any algorithm that does not double-count mergers, 
e.g. a combination of "stitching-3" for short fragmentations 
and "splitting" for long fragmentations, is valid from this per- 
spective. 

The second advantage relates to the timing of the merger, 
and makes us always prefer "splitting" over a combination of 
"splitting" and "stitching". In cases described by the second 
column of Table[T] where both the first, temporary, merger and 
the fragmentation duration are shorter than three snapshots 
and comparable to each other, there's no significant difference 
between "stitching" and "splitting". In such cases it is difficult 
to determine conclusively what the "correct" time is, i.e. that 
of the first or the second merger (but note that Figure |l(a)[ 
as mentioned earlier, suggests that the later merger is usually 
the more physical choice). Nevertheless, when the time spent 
between the fragmentation and the re-merging is larger than 3 
snapshots, the different algorithms perform very differently: 
"stitching-cx)" counts it at the time of the first merger, while 
our method counts it at the time of the re-merger As Figure 
|l(b)| shows, in those cases when the fragmentation phase lasts 
more than three snapshots, it is mostly a few times longer 
than the first, temporary, merger. This further suggests that 
the point of re-merging should be usually considered as the 
more physically appropriate time of merging. 

Our splitting algorithm changes the halo mass function «(> 
M,z) only slightly. In our new halo catalogues there are nat- 
urally fewer massive halos and more small halos, compared 
to the original FOF group catalogues. Figure |2] shows that 
the change does not exceed w 15%, while where the change 
is statistically significant (according to Poisson errors) it does 
not exceeds 3%, for M> 1.5 x 10" M© at allredshifts. This 
difference is smaller than the uncertainty due to e.g. different 
iden tification methods for halos or different simulation codes 
(e.g. lLukic et~ani2007l:lHeitmann et al.ll2007h . 

Our method implicitly assumes that no splits occur after 
z = 0, which is of course physically wrong but technically un- 
avoidable, since the simulation terminates at z = 0. The typi- 
cal time it takes the subhalos to disappear after having merged 
with bigger FOF groups is comparable to the cosmic lookback 
time at z « 0.1. Therefore, starting from z ~ 0.2, it becomes 
impossible to identify a non negligible fraction of the arti- 
ficially linked FOF groups - those that would have split after 
z = had the simulation continued running. In particular, FOF 




Fig. 2. — The ratio between the halo mass function after applying our 
splitting algorithm and the original mass function. The difference does not 
exceed Ri 15%, and where it is statistically significant (according to Poisson 
eiTors, which are shown by the error bars) it does not exceed 3%, for M > 
1.5 X 10" Mq at all redshifts. 



groups at the last snapshot of the simulation are never split by 
our splitting algorithm, since it is impossible to know which 
subhalos of theirs would split "in the future". Our results are 
therefore valid only for z > 0.2, and for use at lower redshifts 
an extrapolation should be used. 

We define the mass of a halo as the mass of all the parti- 
cles gravitationally bound to it, i.e. the sum of its subhalos' 
masses. For w 80%(« 95%) of the halos, the unbound parti- 
cles, which we don't take into account, are less than 5%(20%) 
of the total FOF mass. Also, the halo mass function «(> M,z) 
changes by just a few percent if also the unbound particles are 
included. Nevertheless, including the unbound particles (as 
in e.g. FM08) has a significant effect on the inferred merger 
rate. In ^4.31 we describe this difference, explain its source 
and justify our choice of not including the unbound particles 
in the mass of halos. 

2.3. Definitions of the merger rate: per progenitor halo and 
per descendant halo 

We identify a merger whenever two or more halos at snap- 
shot s have a common descendant at snapshot s+\, and use 
the time/redshift difference between the two snapshots to de- 
fine the merger rate per unit time/redshift, respectively. We 
define a merger as a two-body event, in a way that if n > 2 
halos merge, then n - 1 mergers are recorded, each between 
the most massive one and one of the others. The possibility 
that the mergers occur in a different order, i.e. that some of 
the smaller progenitors merge with one another before merg- 
ing together with the most massive progenitor, is sufficiently 
small so that our results are not strongly affected by it, as 
shown by FM08. 

We derive the merger rate ^tx progenitor halo per unit time 
(or redshift) per mass ratio x: ;y ' , , (x,z,M). Our 

book-keeping is performed as follows. For a merger between 
two halos of masses Mi and M2 we record one merger at mass 
Ml with ratio x = Mi /M2 and one merger at mass M2 with ra- 
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tio X = M2/M1. In this way, all the mergers of each halo with 
any other halo, less or more massive, are recorded. There is 
no double counting in this method, since we are interested in 
the merger rate per halo rather than the absolute number of 
mergers. We attribute a merger to each halo in the pair, and 
each merger event is counted as 2/2 = 1 merger per halo. 

In comparison, for the merger rate per descendant halo, 
only one merger is recorded at mass Mi + M2 with ratio 
X = Ml /Ml (where the indices are defined so that Mi > Mi). 
To better understand the difference, let us consider equal-mass 
mergers (x = 1). The merger rate derived by using the descen- 
dant halo would be one x=\ merger per halo at 2M, while ac- 
cording to our progenitor method it would be one x=\ merger 
per halo at M. Therefore, the typical time scale for a halo of 
mass M to encounter another halo of mass M is the recip- 
rocal of the merger rate per progenitor halo at mass M and 
x=\. The merger rate per descendant halo at mass M can be 
quantitatively very different 03.4l i and has a different physi- 
cal meaning. The merger rate per descendant halo gives the 
time scale on which the population of halos of mass M is cre- 
ated by equal-mass mergers. It is important to keep in mind 
the physical meaning of each of the definitions when imple- 
menting them to physical problems. For example, as we will 
discuss in ^ the appropriate quantity to use when consider- 
ing merger fractions derived via pair counting is the merger 
rate per progenitor halo. 

In the binary merger approximation, which we find to be a 
good one, both definitions of the merger rate are interchange- 
able, given that the halo mass function is known (Appendix 
IaI i. Nevertheless, the merger rate per progenitor halo quan- 
tifies also the rate at which halos merge with more massive 
halos, i.e. the range x < 1, in addition to x > 1. This is im- 
portant because the merger rate per descendant halo is given 
only for the range x > 1 , which represents the rate at which 
halos of mass Mi merge with (or accrete) less massive halos 
of mass Ml /x. However, the range x < 1 represents the rate 
at which halos of mass Mi merge with (or are accreted onto) 
more massive halos of mass Mi /x. To illustrate the impor- 
tance of this difference, consider the way a major merger be- 
tween halos of masses Mi and Mi = Mi /2 is recorded. When 
using the descendant halo, one merger with x = 2 is recorded 
at mass 3 /2Mi = 3M2, while the fact that the halo with Mi 
experienced a major merger is not explicitly accounted for In 
contrast, when using the progenitor halos, two mergers are 
recorded: one with x = 2 at mass Mi and one with x = 1/2 at 
mass Ml. 

3. RESULTS 

3.1. The merger rate per progenitor halo as a function of 
mass ratio 

1 rfWn,, 



TABLE 2 

Fits for the merger rate as a function of mass ratio in Figure 



Figure [3] shows the merger rate 



-{x,z,M) as a 



A'piog-halo dtdx 

function of mass ratio x, for a constant redshift at different 
mass bins (Figure |3(a) and Figure |3(b)| i, as well as for a 
constant mass bin at different redshifts (Figure [3 (c)[ ). In all 
curves, the relation is close to a power law at large x. As x 
decreases, the slope becomes steeper in most curves. At fur- 
ther lower X, in all curves, there is a flattening and the sign of 
the slope changes rapidly as the function starts decreasing to- 
wards very small x. We fit this shape with the following fitting 
function: 
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dN^, 7 A/) = 



Wprog-halo dtdx ' 

AxHl + 1 /xYexpi-ixMjMy') Gyr"' , 



m 



Panel 


Low mass 


High mass 


z 


A 


b 


c 


d 


Mc 




[Mq] 


[Mq] 












[10'" Mq] 


a 


1011.25 


IOII.35 


0.24 


0.006 


-0.18 


0.82 


-0.79 


32290 


a 


IOI2.3 


IOI2.7 


0.24 


0.011 


-0.32 


0.57 


-0.89 


32290 


a 


1013.3 


1013.7 


0.24 


0.012 


-0.23 


0.78 


-1.12 


32290 


b 


1011.25 


IOII.35 


3.06 


0.32 


-0.57 


0.62 


-0.35 


13.4 


b 


1012-3 


10i2.7 


3.06 


1.2 


-0.71 


1.14 


-0.38 


13.4 


b 


10'3 


10"! 


3.06 


2.23 


-0.67 


o<") 


-0.26 


13.4 


c 


IOI2.4 


1012.6 


0.24 


0.011 


-0.32 


0.62 


-0.73 


32290 


c 


IOI2.4 


1012.6 


1.63 


0.1 


-0.51 


0.53 


-0.48 


698 


c 


IOI2.4 


IOI2.6 


3.06 


1.409 


-0.74 


1.37 


-0.39 


13.4 



Note. — (a) For such high mass at high redshift, the 
off affects the merger rate at x large enough that the upturn 
(1 + 1 /xf term is absent. The lowest is actually achieved 
c < 0, but we do not allow that, hence c = is forced. 



exponential cut- 
provided by the 
in this case with 



(1) 



with A = A(z,M), b = b{z,M), c = c{z,M), d = d(z,M) and 
M,=M,(z). 

The x* term describes the shape of the curves for x ^ 1, 
the (1 -I- 1/x)'' term increases the slope at x < 1, and the ex- 
ponential term causes the cut-off at small x, depending on the 
mass of the halos in question and the parameter M^.. We in- 
terpret the exponential cut-off as a consequence of the expo- 
nential cut-off in the mass function: at M/x <C Mc this term's 
contribution is negligible, while at M/x > Mc the function de- 
creases exponentially because there is an exponentially small 
number of halos of mass Mc for halos of mass M to merge 
with. Indeed, there is a close relationship between Mc and the 
dark matter halo m ass at the knee of t he halo mass function 
M, (based on EPS: lMo & Whitell2002l) , with M^ « 30M*. At 
high redshift, the exponential cut-off affects the function al- 
ready at X < 1, so the parameters c, d and Mc become some- 
what degenerate. Therefore, we make use of the relationship 
between Mc and M* to reduce the freedom of the fitting by 
setting log„)(Mf.[ Mq]) « 14.8 - 1 .2z. Table|2]shows the best- 
fitting numerical values for the parameters of the fitting func- 
tion, for the curves that appear in Figure[3] 

To investigate the sensitivity of our results to the time reso- 
lution of the snapshots in the simulation, we have performed 
the following test. For any progenitor snapshot s, we have ex- 
tracted the merger rate skipping the two subsequent snapshots, 
as though the next available snapshot was only s-l-3. At low 
redshift we find a negligible change in the results. At increas- 
ing redshift, Az increases, so when 2 snapshots are skipped, 
there occurs a non-negligible effect of averaging the redshift 
dependence of the merger rate. This is why at z « 3 a decrease 
of up to w 20% can be seen for this "skip 2" method, a differ- 
ence that is consistent with being just the result of averaging 
the redshift dependence. This Az-convergence is comparable 
to the one found by FM08. 

As can be seen in Figure [3] and Table |2] the shape of 
_i — '^'^"j'i" (x 7 M) changes with mass and redshift, so that 

A'prog-halo utax 

no global values for the fitting parameters A, b, c & d can 
be obtained. This reflects deviation from self similarity, in 
the sense that halos of different masses and at different red- 
shifts have different weights for merging with halos of differ- 
ent mass ratios. At a given redshift (Figure [3(a)| and Figure 
|3(b)| i, more massive halos have slightly more minor mergers 
(x > 1, consistently with FM08) while they are being accreted 
onto larger halos (x <C 1) at a much lower rate than less mas- 
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log^jj(M[M^J)=11.3 
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Fig. 3. — The merger rate 



1 



;rgcr 



, , -{x,z,M) as a function of mass ratio x. Panel (a) displays the rates for several progenitor masses, at fixed redshift 

z ~ 0.24, and Panel (b) shows the same for z~ 3. Panel (c) displays the rates at different redshifts, for 10'^ Mq < M < 2 x lO'- Mq halos. Asterisks show 
data based on the simulation, and the solid curves show our fit based on equation (T) and Table|2] To guide the eye, two veilical lines show the range of x within 
which we define major mergers. We note that x is limited at large values by the mass resolution limit of the simulation: we show only mergers between M and 
M/x > 2 X M,„„, = 4.72 X lO'" Mq. At small values x is Hmited by the simulation box size. The bins in x are logarithmically equally-spaced, except for at the 
left-most side of each curve, where bins were constructed so that they contain no less than 5 mergers. 



sive halos. Similarly, for a given halo mass (Figure [3(c)| l, at 
high redshift the minor merger rate (per unit time) is higher, 
while the rate of being accreted onto bigger halos is lower. 

3.2. The major merger rate per progenitor halo as a function 
of mass and redshift 

For a quantitative comparison between similar masses at 
different redshifts or different masses at fixed redshifts, a spe- 
cific value or range in x must be chosen. We investigate the 
dependence on mass and redshift of the possibly most inter- 
esting range: 1/3 < x < 3, i.e. the major merger rate: 



that equation (O is valid for 

M < M,(z)/5 w 6M,(z) « 10^4-^-'-2' Mq. (4) 

Similarly, for a given mass M, at low enough redshifts 
the major merger rate increases with increasing redshift, 
as shown in Figure |4(c)| The relation is of the form 
^(z,M) cx (1 with /? « 2. But since M,(z) 



1 



dN, 



major— merger 



Nr 



prog-halo 



dt 



1 



dN„ 



1 /3 -^pi-og-halo 



(x,z,M)ci(e) 



dtdx 

Note that, following the discussion in ^2.31 in cases where in- 
tegrating over the range 1 /3 < x < 3 counts the same merger 
twice, also both halos are counted in A'prog-haio- In other 
words, each merger always contributes exactly one count to 
the merger rate per halo, so there is no double counting in 
equation (|2|i. 

Figure |4(a)| shows the dependence of 

1 rfA-majcr-magcr^^^^-) ^^^^ at Z « 0.24. At 



dt 



^prog-hiilo 

M < 10'^'^ M0 the major merger rate increases with 
increasing mass, and a relation of the form 



1 



dN, 



major— merger 



dt 



(z,M) cx log(M) + a 



(3) 



*prog-halo 

holds, as the solid green line in Figure |4(a)| shows. How- 
ever, for higher masses the major merger rate decreases 
with increasing mass, because the exponential cut-off of 
_i — dN,„,,a y-. gjj^^jj ^ affects even the range x > 1/3. 

A'prog-haio dtdx b I 

At each redshift the critical mass above which equation (O 
no longer holds is different. At high redshift there are less 
massive halos, so the exponential cut-off occurs at x = 1/3 al- 
ready for lower masses, compared with low redshift. We find 



dN^: 



A'prog-haio dl 

decreases with increasing redshift, there is a redshift above 
which M > 6M,(z). At that redshift the major merger rate be- 
gins to be affected by the exponential cut-off, and the power 
law approximation breaks down. The major merger rate starts 
decreasing with increasing redshift only at much higher red- 
shifts (e.g. while 6M,(z « 0.9) = 10'^ Mq, the major merger 
rate of 10^^ halos starts decreasing with increasing red- 
shift only at z > 3). We note that the regime that is affected by 
the cut-off exists for the merger rate per progenitor halo that 
we present in this Section. It does not appear in the merger 
rate per descendant halo, as will be shown in 33.3l and ^3.41 

As long as equation (|4| holds, we find the following fitting 
function for the major merger rate as a function of mass and 
redshift: 



1 



major— merger 



prog-halo 



dt 



(z,M)=/?(^io + a)(l+z)^, (5) 



with /? = 0.0075 Gyr"', a = 1 ±0.1, /3 = 2.1 ±0.1 and 
logi()( 



M 



lOV i()io„ 



). It is accurate to within « 10%, an inaccuracy 

that is included in the errors quoted for the fitting parameters. 
Because of the inability to properly split halos in the very last 
snapshots of the simulation at z < 0.2 (as explained in ^2.2l i 
we do not include this range in our fitting. 

An illustration of the regime where equation (HI does not 
hold is shown in Figure [3tb)| and Figure [4(b)| Since M^z = 
3) ^ lO"'^ M0, all the mass ranges we can probe at this red- 
shift are in the regime where equation (|4| doesn't hold. In 
Figure |3(b)| it is shown that at z ~ 3 the exponential cut-off 
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Fig. 4. — The major merger rate per progenitor halo as a function of mass at z~ 0.24 (a) and at z Ri 3 (b), and as a function of redshift for 10 Mq < M < 
IO'^^Mq (c). The asterisks show our analysis of the simulation, while the solid green curves show the fitting function equation (5). The fitting functions fit 
the data well for low enough values of mass and redshift so that equation (4) holds. In this regime the major merger rate per unit time increases steeply with 
increasing redshift and mildly with increasing mass. Above the threshold set by equation j4), we find the mass/redshift dependencies to break. Specifically, at 
high redshift (b) the whole mass range available from the simulation is above the exponential cut-off mass given by equation j4), therefore the major merger rate 
is seen to almost always decrease with increasing mass. Equation (5) is not vaHd in that regime. 




iog,„(M[M^„j) 



Fig. 5. — The major merger rate per descendant halo as a function of mass atz~ 0.24 (a) and at z Ri 3 (b), and as a function of redshift for 10 Mq < M < 
10'2 6 (c) xhe asterisks show our analysis of the simulation and the red curves show the merger rate as quantified by FMOS's fitting formula (and integrated 
over the major merger range of mass ratios, between x = I and x = 3). A comparison to Figure |4] shows that even above the threshold set by equation j4) the 
mass/redshift dependencies of the merger rate per descendant halo do not change, in contrast with the merger rate per progenitor halo. A comparison to FMOS's 
results shows that the major merger rate is lower once our splitting algorithm and our mass definition are used. 



occurs at much larger x (compared with Figure [3(a)]), so that 
also the major and minor merger regimes are affected. This 
means that the merger rate is lower for more massive halos 
even at > 1 , as opposed to low redshift, where this is true 
only for x < 0. 1 . The major merger rate, explicitly as a func- 
tion of mass, is shown in Figure [4(b)| 

It is interesting to note that the major merger rate per unit 
redshift, 

1 dN major— merser / . .^^ 1 ^/A^ major— merger , . 

-1 —iz,M) = — — ^(z, A 

^*prog-halo ^*prog— halo 

is approximately constant with redshift for the range that is 
not affected by the exponential cut-off, since j;iz) approxi- 
mately cancels out the redshift dependence in equation (|5]l. 
At redshifts higher than the break redshift, the major merger 
rate per unit redshift decreases steeply with increasing red- 
shift. The bimodal fit: 

0.38 ±0.02,z< 2.5, 



0.38-0.063(z-2.5),z>2.5 (7) 

describes the major merger rate per unit redshift for 10'^ Mq 
halos as a function of redshift. It can be applied for different 
masses using equation (|5]l to scale with mass, and equation 
(in to find the "break" redshift. 

If the mass dependence of the major merger rate, which is 
weak, is neglected, equation ^ can be integrated between z, 
and Zf < Zi to achieve the average number of major mergers 
throughout the formation history of halos. If a population of 
halos of mass M is chosen at redshift z, where M < 6M*(z,), 
then for z < Zi the major merger rate per unit redshift is 
roughly constant with redshift. Therefore the very simple lin- 
ear relation for the average number of major mergers halos of 
mass M at redshift z, will undergo until redshift z/ is 

.per-prog iZi,Zf,M) ~ 0.13 X (pio+ l)(z,-z/), (8) 
while the number of major mergers individual halos undergo 
is of course distributed around this average. 
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3.3. The merger rate per descendant halo 

We find that the merger rate per descendant halo can be fit 
much better with a global fitting formula, compared with the 
rate per progenitor halo. We adopt the fitting form of FM08 
(albeit keeping our mass ratio variable x = 1 /^), and give dif- 
ferent best-fitting parameters that express our different treat- 
ment of fragmentations i ^2.2\ and different mass definition 
(El: 



1 



dSr 



^I^(x,z,M) =AM'^2^''exp((x/xry-^, (9) 

A'desc-haio dzdx dz 

whereMi2 =M/1O'^M0 and Sc « 1. 69 /D (z), which we esti- 
mate using the approximation provided by iNeistein & Dekell 
(|2008a). 

Our best-fitting parameters are: A = 0.06, a = 0.12, b = -0.2, 
X = 2.5 and 7 = 0.5. Note that with the mass ratio definition 
used by FM08 our parameters correspond, in their notation, 
to: A = 0.06, a = 0.12, /3 = -Z7-2 = -1.8, |= 1/f =0.4, 7 = 
0.5, 77 = 1 and M = 10'^ Mq. A detailed comparison with the 
results of FM08 is presented in ^4.11 

We find this formula to fit the merger rate per descendant 
halo per unit redshift with deviations of up to « 20% for all 
the mass range probed by the Millennium Simulation at red- 
shifts z < 4. One systematic exception is at 0.5 < z < 1.5 
and 30 < X < 1000, where the fitting formula tends to over- 
estimate the merger rate by up to 50%. At z > 4 the redshift 
dependence as well as the mass dependence become stronger, 
and we do not make an attempt to fit that regime. 

Also the merger rate per descendant halo can be integrated 
across cosmic times (neglecting its mass dependence) to ob- 
tain the mean number of major mergers that halos of mass 
M at redshift Zf have undergone since redshift z,. If we ap- 
proximate ^ as 1.25 at z > 1 and 0.8 + 0.32z at z < 1 (an 
approximation that holds with deviations < 10%) then we ob- 
tain 

A'MM,per-desc(Zi,Z/,M) « 0.43MJ'2 X 

'(Zi-Zf) if4>z,-,z/>l 
z,-0.65z/-0.14z}-0.21if 4 > zi > l,z/ <1D) 
,0.65(z,-z/ + (z2-z})/5) ifz,,z/< 1. 

3.4. Comparison between the merger rate per progenitor 
halo and per descendant halo 

For mergers of x « 1 at high mass/redshift the different def- 
inition of the merger rate, i.e. per progenitor halo versus per 
descendant halo 02.3b . is an important issue. Where equa- 
tion (01 is no longer satisfied, the functional behaviour of the 
two is different: defining the merger rate per progenitor halo, 
the exponential cut-off of the merger rate enters into the x > 1 
range, whereas defining it per descendant halo, the function 
increases all the way to x = 1, beyond which x is no longer de- 
fined. This difference is caused by attributing major mergers 
(e.g. between identical halos of mass M) to their approximate 
sum of masses. The same number of mergers is recorded by 
both methods, but there are far fewer halos with mass 2M than 
halos with mass M in the exponential cut-off regime. 

For X 3> 1, or the minor merger regime, the difference be- 
tween the two definitions is small. For example, the integrated 
rate of mergers with 10 < x < 100 per progenitor halo is lower 
by « 10% at z < 1 and M « lO'-' M© than the same rate per 
descendant halo. 

We further consider the major merger rate, i.e. the merger 
rate integrated over l/3<x<3(orl <x<3 for the merger 



rate per descendant halo). Figure |5] shows the merger rate 
per descendant halo, keeping our halo mass definition and 
treatment of halo fragmentation {blue, asterisks). Comparing 
Figure |5] to Figure |4] illustrates that in the low mass/redshift 
regime where equation (|4]l holds, the major merger rate is only 
slightly affected by the merger rate definition. As described 
in ^3.21 above the threshold set by equation (|4|, we find the 
mass/redshift dependencies of the merger rate per progenitor 
halo to break, a different behaviour than that of the m erg er rat e 
per descendant halo. At high redshift (Figures [5(b)| and r4(b)[ ), 
where all the masses that can be probed using the Millennium 
Simulation are above the "break mass" (equation (|4|), the dif- 
ferent definitions lead to very different results at almost all 
masses. 

Finally, only the merger rate per progenitor halo quantifies 
the rate of "destruction" events, where halos are accreted onto 
more massive halos, in an explicit way. This "destruction" 
merger rate, x < 1, is not defined by the merger rate per de- 
scendant halo alone. To obtain the merger rate with x < 1 from 
the merger rate per descendant halo, the halo mass function 
must be known. In Appendix lAl we develop analytic formu- 
lae that convert merger rates per descendant halo into merger 
rates per progenitor halo and vice versa, and also relate the 
merger rate per progenitor halo at x < 1 to that at x > 1 . The 
conversion formula from the merger rate per descendant halo 
to the merger rate per progenitor halo is given by 



Rp(Zp,M,x > 1) = 
Rdizd, -^M^xy^Nhizd, '^M)/Nh(zp,M), 



(11) 



RpiZp,M,x< 1) = 

RAzdA-,+ m, + mh(zd,(j + m)x-yNh(zp,Mj,i2) 

where Rp is the merger rate per progenitor halo, R^ is the 
merger rate per descendant halo, Ni, is the mass function de- 
fined as Nh(M) = , and Zp and Zd are the progenitor and 
descendant redshifts, respectively. We verified numerically 
that those formulae describe the relations between the differ- 
ent merger rate definitions well. 

4. COMPARISON TO PREVIOUS WORK 
4.1. Comparison to FM08 

In Section 12.21 we have described the differences between 
our method and FM08's method for building new halo trees 
and defining the halo mass. We show in this section that those 
differences have a significant effect on the results. We first 
discuss the differences between the dependencies on mass ra- 
tio, then compare integrated minor and major merger rates, 
and finally indicate the different time and mass evolutions. 

Figure |6] shows the merger rate per descendant halo as a 
function of mass ratio for 10'"* Mq halos at z ~ 0.4. For mi- 
nor mergers FM08* also find that the merger rate varies as a 
power law of the mass ratio. However, they find the power 
law index to be « (/3 = -2.01), i.e. find more minor merg- 
ers. This can be seen at x ^ 1 in Figure |6] The difference 
arises because FM08 accept most temporary links between 
FOF groups as mergers. This leads to an artificial inflation of 

* FM09's "splitting-3" method gives an almost identical merger rate to 
FMOS's fiducial "stitching-3" method, in spite of the differences between the 
algorithms. An examination of Table [T] shows that it is probable that those 
differences compensate each other between the scenarios in the third and fifth 
columns of that Table. 
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Fig. 6. — The merger rate per descendant halo for 10 Mq halos aX z~ 
0.4, comparing the results of our splitting algorithm to FMOS's "stitching" 
algorithm andFM09's "splitting-3" algorithm (C.-P. Ma, priv. comm.), which 
both give virtually the same results. We find a lower merger rate, starting from 
f» 6% at X < 10 and increasing steadily to a factor of 2 at ,t 1000. Abou t 
half of the difference originates in the different halo mass definition ( j|4.3t , 
and about half from our algorithm that never counts the same merger more 
than once. 
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Fig. 7. — The integrated merger rate per descendant halo per unit time of 
minor mergers between 10:1 and 100 : 1, as a function of mass at z « 0.24. 
The mass resolution of the simulation allows this rate to be directly mea- 
sured only forM > IO'^ '^Mq. Appropriate treatment of halo fragmentation 
and an appropriate halo mass definition are significant in this regime. Not 
giving any special treatment to the fragmentation of halos, thereby counting 
artificial mergers (FMOS's "snipping" method, dashed red) leads to a signif- 
icant overestimate of the merger rate. FMOS's "stitching" method (solid red) 
still overestimates this minor merger rate by a factor 1.3 compared to our 
method {green, asterisks), due to their acceptance as mergers of most tempo- 
rary links between FOF groups and their definition of a halo mass. 



the number of minor mergers. This was anticipated by FM08 
themselves, and is also supported by the finding that subhalos 
that a re ejected out of their host halos are preferentially of low 
mass dLudlow et alJ l20Q8). The second significant parameter 
that depends on the tree building method is x (see equation 
Q). For the "snipping" method FM08 found that the best- 



fitting value is w 58, and for the "stitching" method « 10, 
while we find the best value to be x = 2.5. This means that the 
dependence of the merger rate on mass ratio in our method is 
closer to a pure power-law (since the exponential term affects 
mostly X < x). The larger the number of false (mostly minor) 
mergers that are counted, the stronger the exponential term 
should be, and the steeper the power-law. 

The integrated minor merger rate per descendant halo of 
mergers between 10 : 1 and 100 : 1 at z « 0.24 as a function 
of mass is shown in Figure |2l It shows that the "snipping" 
and "stitching" methods {dashed and solid curves) overesti- 
mate the minor merger rate that we find. Roughly half of 
this difference originates from properly removing all artificial 
mergers, and equally important is the difference in the mass 
definition of halos (see 34.3l l. Note that while Table [T] shows 
that w 80% of all mergers are treated equally by all methods, a 
large fraction of all mergers are rather major mergers between 
rather low mass halos, simply because of the large abundance 
of low mass halos. Among minor mergers, which involve ha- 
los much more massive than the mass resolution limit, many 
belong to the other w 20% that are treated differently by the 
different methods. 

The comparison in the major merger regime can be seen 
in Figure |5] The merger rate we find is systematically lower 
than FM08's also due to our splitting algorithm but, in this 
regime, mainly due to the different halo mass definition. This 
is especially true at low mass (Figure [5(a)[ ) and low redshift 
(Figure [5(c)l i. 

We find the redshift dependence of the merger rate to be 
proportional to ^ , while FM08 find a slightly weaker red- 
shift evolution at low redshift by introducing (^)'' with 
?7 w 0.3. The consequence of ?/ = 1 is discussed further in 34.21 
We also find a slightly stronger mass dependence, a = 0.12 
rather than a « 0.08. 

4.2. Comparison to the EPS model 

Figure [U compares different results for the merger rate per 
descendant halo of 10'"* halos at z « 0.4. The merger rate 
per descendant halo is the natural quantity that is obtained 
from the EPS model, and does not require explicitly the mass 
function of halos. 

FM08 found that the EPS model predicts more major merg- 
ers and (especially) fewer minor mergers compared with the 
Millennium Simulation. The slope of the merger rate at x ^ 1 
is determined by the parameter b in equation (|9]l. The larger 
it is, the more minor mergers dominate the number of merg- 
ers. Indeed, FM08's "snipping" method, which doesn't re- 
ject any artificial mergers, finds b = 0.17, while their fiducial 
"stitching" method (red curve in Figure [8(a)] i finds = 0.01. 
The value we find, by rejecti ng all artificial mergers (blue 
curve with asterisks in Figure 8(a)| i, \s b = -0.2, while the 
iLacev & Cold (|1993!) EPS prediction (green curve with open 
circles in Figure [8(a)] l is « -0.5. Therefore, we find, like 
FM08, that the EPS model overpredicts the major merger rate, 
but we find a better agreement with EPS in the minor merger 
regime, because we don't acc ept artificial minor mergers. 

Recently Neistein & Dekell ( l2008b) have constructed a new 
method for predicting merger rates from the EPS model, by 
avoiding the assumption of binary mergers, which leads to 
inconsistency withi n EPS. They find, in agreement with the 
simulation, that the lLacey & Cold (11993) merger rates are too 
low in the minor merger regime. Moreover, we find that the 
iNeistein & Dekell (l2008bh merger rate (green curve with filled 
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circles in Figure [8(a)] i has atx 1 a slope of /? = -0.19 ±0.02, 
i.e. its shape agrees remarkably well with what we find in 
the Millennium Simulation. The iNeistein & Dekell (l2008bh 
merger rate is higher than the merger rate we find in the Mil- 
lennium Simulation by a factor of w 1.5 ±0.3, for all mass 
ratios, halo masses and redshifts, with almost no systematic 
dependence on any of those parameters. 

This result is further d emonstrated in Figure |8(b)| The 
INeistein & Dekell (l2008bl) merger rate is constant with respect 
to the natural dimensionless EPS time variable Sc, which is 
exactly the same dependence we find in the simulation. Fig- 
ure [8(b)] compares the major merger rate per descendant halo 
of M « 10'^-^ M0 halos between FM08's "stitching" method, 
the method presented here, and INeistein & Dekell (l2008bl) . 
The redshift dependence of FM08's method is significantly 
different at low redshift (originating numerically from their 
best-fitting value of 77 « 0.3), because artificial mergers ap- 
pear preferentially at low redshift. The redshift dependence 
of the merger rate in the Millennium Simulati on based on our 
method matches well that of EPS (see also Neist ein & Dekell 
(|2008a)), but the normalisation, as already mentioned, has an 
offset. 

4.3. The role of the halo mass definition 

There is a substantial uncertainty as to how to determine 
the boundaries of halos in N-body simulations. Sometimes 
the mass of halos is taken as the mass inside a sphere, within 
which the density equals the expected density of virialised 
groups in the spherical collapse model. This class of def- 
initions is not well-suited for our purposes, because they 
are less reliable f or halos undergoing mergers (White 200 It 
iLukic et al]l2008l) . Particles grouped together by the EOF al- 
gorithm are considered to correspond to virialised dark mat- 
ter halos, because their average density approximately equals 
the expected density of virialised groups in the spherical col- 
lapse model, yet the halos can have any shape and are not 
assumed to be spherical. On the other hand, undesir ed ef fects 
like particle bridges between halos (as discussed in ^2.21 ) and 
spurious linkage of particles to groups are inherent to the algo- 
rithm. An advantage of SUBFIND over EOF is that it subjects 
the (geometrically-identified) groups of particles to a dynam- 
ical test. Only particles that are found to be gravitationally 
bound to their subhalo are included as part of their subhalo's 
mass. Therefore, two reasonable definitions we consider for 
the mass of a halo are: (1) the total mass of all the particles as- 
sociated with the halo by the EOF algorithm (as in e.g. FM08), 
and (2) the total mass of all the particles that belong to the halo 
and are also gravitationally bound to any of its subhalos. 

As described in ^2.21 for the population of halos as a whole, 
the difference between those two definitions is not large. For 
example, the average fraction of unbound particles is « 3%. 
Nevertheless, there is a distinct population of halos for which 
the typical fraction is much larger These are halos that are 
about to undergo a significant merger, i.e. undergo a major 
merger or be accreted onto a more massive halo. First we 
describe how this affects the merger rate, and then interpret 
this phenomenon and justify our choice of not including the 
unbound particles when computing the merger rate. 

Halos at snapshot s that are about to undergo a merger have 
on average higher ratios of total FOE mass over bound SUB- 
FIND mass. We find that the total-over-bound mass ratio of 
a halo correlates with two quantities related to a halo's next 
merger: 1) the time before the next merger begins, relatively 
to the dynamical time oc H{z)~^, and 2) the mass ratio of that 



merger. Specifically, as a halo approaches a merger with an- 
other halo of comparable or larger mass, its total-over-bound 
mass ratio increases with time, an effect that is stronger as the 
mass of the other halo is larger. This phenomenon can be well 
quantified for x <= 3 by 

Mfof 



-(x<= 3)=1.02xf, 



-0.018 
leff.dyn 



X X 



-0.005 X 



Ieft.dyn 
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Mbound 

where ti^ft^yn = tieftH(z), and ti^f, is the time left before the 
merger starts (which is defined for this purpose as the mid- 
dle of the time interval between the adjacent snapshots of the 
progenitors and the descendant). denotes the geomet- 

'''hoiiiul 

rical mean over all halos, a quantity that is used because the 
scatter is large. Halos that are approaching a minor merger 
(jc > 3) show almost no enhancement of the total-over-bound 
mass ratio, therefore we fit it as 



M, 



™^x>3)=1.02xr^«lJJ„. 



)l/f ^ ' it'ji.uyn- (14) 

Mbound 

Figure|9]demonstrates this phenomenon and its description by 
equations ( fT3] l and ( fT4l l for halos at « 3. 

Since the more massive halo of the merging pair (or group) 
retains the typical value of sa 1.03, and the less massive ha- 
los have untypically large total-over-bound mass ratios, the 
mass ratios of mergers shift towards smaller values (merg- 
ers become more equal-mass) when the total mass is taken 
into account, compared with the choice of the bound mass as 
the halo mass. The effect on the merger rate is pronounced. 
When the halo mass is taken as the total mass, there is a de- 
ficiency of very high mass ratio mergers near the resolution 
limit, but since we restrict our analysis to mergers with ha- 
los of M/x > 2 X M,ni„, this is not seen. What is seen is an 
enhancement of the merger rate for every jc < M/(2 x M„„„). 
The merger rate using the total EOF mass is larger by typi- 
cally 20% and up to w 50% than our fiducial method. Fig- 
ure [TO] demonstrates this difference. Furthermore, Figure [TOl 
shows that using the mean relations equations ( fTJi l and ( fT4l i 
to interchange between the two halo mass definitions allows 
reproducing the different resulting merger rates. 

The total-over-bound mass ratio of halos correlates also 
with the environment in which they reside. Figure [TTI shows 
that halos in denser environments have higher total-over- 
bound mass ratios {solid curves). This trend holds also among 
halos that are about to undergo a x < 1 merger in the following 
snapshot {thin blue). Yet, once each halo's total-over-bound 
mass ratio is normalised by the value expected for it by equa- 
tions ( fTST l and ( fT4b . the correlation with the environment al- 
most disappears {dashed curves). The difference, at a given 
overdensity, between all halos and halos about to undergo a 
x< \ merger shows that it is not possible to use the correlation 
with the environmental overdensity to disentangle the depen- 
dence shown in Figure |9l This means that the environmental 
dependence is probably just a second-order correlation, orig- 
inating from the correlation shown in Figure |9] in combina- 
tion with the correlation between mergers and environment 
(cf. FM09). 

The strong correlation of the total-over-bound mass ratio of 
halos with their proximity to their next merger and its mass ra- 
tio, and the fact that there is almost no additional dependence 
on environment, suggest that the mergers themselves are re- 
sponsible for the change in the total-over-bound mass ratio. 
The question to be asked is whether this is a gravitational 
effect, which changes the bound mass even before the FOF 
groups merge, or a numerical effect, which changes the FOF 



The Millennium Simulation halo merger rate 



11 




Fig. 8. — A comparison between different results for tlie merger rate per descendant lialo per unit time. Panel (a) illustrates the merger rate for M K, 10''' Mq 
at z 0.4. Four different methods are shown: using the splitting algorithm and the mass definition presented in ^(hlue, asterisks), FMOS's fit to their "stitching" 
method (dashed red), the Neistein & Dekel (2008b) EPS method (green, filled circles) and the Lacev & Cole 1 1993) EPS method (dashed green, open circles). 
At X < 30, the two EPS methods are similar and exceed the two methods based on the Millennium Simulation by 70%. As x increases towards more minor 
mergers, FM08's merger rate is the shallowest, finding the highest number of minor mergers, while the EPS method of Lacey & Cole 1 1993) is the steepest, 
finding the lowest number of minor mergers. On the other hand, the slopes found by o ur analysis and by the EPS method of Neistein & Dekel (2008b) are 
similar, so that a roughly constant ratio is kept between the inferred merger rates ( j|4.2t . Panel (b) compares the integrated major merger rate (1 < .x: < 3) of 
M lO'^^ Mq halos, as a function of redshift, that results from our method (blue, asterisks), FM08's "stitching" method (dashed red) and the Neistein & Deke| 
(2p08b) EPS method, with which both other methods are normalised. The redshift dependence resulting from our analysis of the simulation agrees well with that 
of the EPS prediction (yet, again, with an offset in the normalisation), while FM08 find more mergers at low redshift because their method is more sensitive to 
artificial mergers, which are more common at low redshift. 
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Fig . 9 . — The ratio of the total EOF mass of halos to their total SUBFIND 
bound mass, as a function of the mass ratio of the merger they are about to 
undergo. The more massive the other halo (the smaller .v), the larger the total- 
over-bound mass ratio enhancement, but minor mergers (x > 3, to the right 
of the horizontal line) have almost no effect. The more imminent the merger, 
normalised to the dynamical time, the larger its influence on the halo's total- 
over-bound mass ratio. These trends are shown by 4 different curves, each 
representing halos that are about to undergo such a merger in a fixed "nor- 
malised time" into their future, as indicated by the legend. All halos shown 
here are at the same snapshot, z ~ 3, but the fitting formulas equations U3) 
and )14t , which are shown by the dashed curves, hold for all redshifts. Since 
H(z f» 3) ~ 0.3 Gyr"', the curves in this figure coirespond to halos that are 
about to undergo a merger in 0.1,0.3, 1.2.4 Gyr into their future, which cor- 
responds to 1,2,5, 10 simulation snapshots, from top to bottom respectively. 



group mass even before the gravitational interaction is signifi- 
cant. To give a fully satisfactory answer to that question, a de- 
tailed dynamical analysis is needed, which cannot be obtained 
with just the merger trees, and is outside the scope of this pa- 
per Nevertheless, Figure [12] shows that it is not the bound 
mass that decreases before the merger, but the total mass that 
increases more rapidly, that causes the total-over-bound mass 
ratio to increase. We postulate that a significant amount of 
particles are added to the outskirts of the FOF group during 
the few snapshots before the merger, but are not found to be 
bound to it. Therefore, we favour the interpretation of the 
bound mass as the "true" mass of the halo. 

5. COMPARISON TO OBSERVATIONS 

The two different definitions of the merger rate, per pro- 
genitor halo and per descendant halo, have different physi- 
cal meanings, and they correspond to the two different ob- 
servational approaches towards measuring the galaxy merger 
rate. When the merger rate is measured via pair counting, 
the mass (or luminosity) of each galaxy in the pair is mea- 
sured separately, and the merger fraction/rate is attributed 
to the measured progenitor galaxy mass. The merger rate 
in this case is related to the time scale for galaxies of a 
given mass/luminosity to encounter other galaxies of a given 
mass/luminosity, analogously to the merger rate per progen- 
itor halo. In the case of morphologic al/kinematical identifi- 
cation of disturbed galaxies that show signs of mergers, it is 
difficult to infer the mass of each of the original components 
that have merged, and the merger is then attributed to the total 
mass of the system, i.e. the descendant mass. In that case, the 
merger rate is related to the time scale on which a population 
of galaxies is created by mergers. While such observations 
still suffer from large uncertainties, it is important to realise 
that the merger rates that are inferred by both methods are not 
the same quantity. Those two quantities are related to each 
other by our conversion formula, in a fully analogous way to 
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Fig. 10. — The merger rate per progenitor halo for halos of M 10'^ Mq at 
z~3, extracted using the total FOF group mass as the halo mass {solid green) 
and using only the bound mass as the halo mass (dashed red). When the halo 
mass is defined as the total FOF mass, including gravitationally unbound 
particles, the inferred merger rate is higher by up to 50% (in a similar way 
for all masses at all redshifts). This is caused by the fact that halos that are 
about to merge with more massive halos have unusually high total-to-bound 
mass ratios. To check whether this is the only source of difference, we used 
the merger trees in which th e ma ss is defi ned as the total mass, and changed 
them according to equations il3i and U4K computing for z~i and a 

1 snapshot difference. Specifically, each halo's mass was reduced by 4.55%, 
except for the mass of halos about to undergo a merger of mass ratio x <= 3, 
which was reduced by a factor of 1 .OSTx"^"^*** . Although this model doesn't 
include the scatter in the total-over-bound mass ratio, the result (dotted blue) 
is very similar to the results achieved directly with each halo's bound mass. 
The integral of the merger rates with both mass definitions between x = and 
X = M/Mmin is by construction equal. Nevertheless, we restrict our results to 
X < M/(2 X M„ii„), so the sharp, artificial, drop of the solid green curve at 
X > M/(2 X Mmin) is not seen. 



the two definitions of the merger rate of dark matter halos. 

Some evidence for the expected difference between the 
two method s can be found in the literature. For example, 
iMaller et al.1 (12006,) measure the galaxy major merger rate in 
a cosmological hydrodynamic simulation, and define it as be- 
ing per descendant galaxy. They find a steep monotonic re- 
lation between increasing mass and increasing merger rate, 
simila r to the trend of the h alo merger rate per descendant. 
Also Conselice et al.l (l200 8i). measuring the merger fraction 
using a morphological investigation, find that the merger frac- 
tion (at z > 1.5) increases strongly with increasing mass. On 
the other hand, Patton & Atfield (2008) find that the major 
merger rate of galaxies, obtained via observed pair counting, 
peaks with respect to luminosity, therefore also with respect 
to mass. This different dependence, which at first might seem 
to be at odds with th e findings of 'Mailer et al.l (l2006l ) and 
IConselice et al.l ( l2008h . is ac tually qu alitativel y expe cted once 
it is taken into account that iPatton & Atfieldl (l2008l) measure 
(implicitly) the merger rate per progenitor galaxy. Therefore, 
when comparing either observed or simulated galaxy merger 
rates, the nature of the observations and the definition used to 
analyse the simu lations must be taken into account. Similarly, 
IConselice et al.l (*2008) find that the redshift dependence of the 
merge r fraction becomes st ronger for higher mass galaxies, 
while Ide Ravel et al.l (l2008l) . measuring the pair fraction of 
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Fig. 1 1 . — The total-over-bound mass ratio of FOF groups as a function 
of their environment and merging state. Halos residing in denser environ- 
ments have on average higher total-over-bound mass ratios (thick red). The 
same is true for a sub-population of halos that are just about to merge onto 
a more massive halo (thin bhte). Also, for a given overdensity, halos that 
are just about to merge onto a more massive halo have significantly higher 
total-over-bound mass ratios, in accordance with Figure Yet, once each 
halo's total-over -bou nd mass ratio is normalised by the value expected for 
it by equations )13t and )14t , the correlation with the environment almost 
disappears (dashed curx'es). This figure shows halos at z « 1, but we find 
the same trends at any redshift, and for any available choice of a smoothing 
length: 1.25,2.5,5, 10 Mpc. 



galaxies, find exactly the opposite trend, i.e. that higher mass 
galaxies have a shalower redshift evolution. This is again in 
qualitative agreement with our expectations based on the dif- 
ference between merger rates per progenitor and per descen- 
dant halo/galaxy. 

To quantitatively compare observations with the frequency 
of galaxy mergers predicted by simulations, a treatment 
of baryonic physics must be included. Such a compari- 
son is outside the scope of this p aper, but can be found 
in e.g. IBertone & Conselic3 (12009 *) and P. F. Hopkins et 
al. (2009, in preparation). However, we make here a compari- 
son between the halo merger fraction (per descendant halo) 
in the Millennium Simulation and galaxy merger fractions 
(per descendant galaxy) from observations. In Figure [13] we 
show a compilation of observed galaxy merger fractions as a 
function of redshift. Two features are apparent: the merger 
fraction increases with redshift, and the scatter at any given 
redshift is roughly a factor of 5. If we distinguish between 
fractions according to a rough luminosity/mass criteria, the 
large scatter remains. This is because the scatter is caused 
by numerous factors, e.g. selection of different populations 
and different wavebands, different techniques for identifying 
mergers, as well as cosmic variance. 

To extract merger fractions from the simulat ion, a time 
scale must be associated with mergers (see also ' Genel et aP 
(2008)). We use the dynamic al friction time scale f o r dark 
matter halo mergers found by iBovlan-Kolchin et al.1 (l2008h 
based on merger simulations. We average the orbital param- 
eters and thereby obtain Tmerger = 0-'^ fa(i+)-) ^hx-^'' ^y^^ where r 
is the mass ratio and H(z) is the Hubble constant at redshift 
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Fig. 12. — Mean mass growth histories for halos that are about to un- 
dergo a merger of 0.0005 < A" < 0.01 at z Si (dotted) or z » 1.5 (solid), for 
10 subsequent snapshots before the merger. Both the total FOF group mass 
(empty circles) and the bound mass (filled circles) are shown. Each curve is 
normalised to 1 at the time of 10 snapshots before the merger. It is seen that 
the total mass grows faster than the bound mass, which results in increased 
total-over-bound mass ratios as these halos approach a merger with a more 
massive halo. For halos not about to undergo a significant merger, the total 
and bound mass grow at the same rate, thereby conserving the mean ratio of 
^ 1.03. 



z. Whenever a major merger occurs between snapshots s and 
i + 1 02.3I I. we tag the most massive progenitor at snapshot s, 
and its descendants in following snapshots, as "undergoing a 
major merger" for a time T„,erger after snapshot s. This allows 
us to determine, for any snapshot and halo mass range, the 
fraction of halos that are instantaneously undergoing major 
mergers, i.e. the major merger fraction. In Figure[T3]we plot 
the resulting major merger fraction as a function of redshift 
for three halo masses. 

We note that by tagging the progenitor halo only at snap- 
shot s and the descendant halos in the following snapshots for 
a duration T„erger, we derive a merger fraction that is mostly 
"per descendant halo". We do that for two reasons. First, 
while a halo that is a descendant of a merger is a halo that 
hosts galaxies that are about to merge, it is probably obser- 
vationally less relevant when two halos are still approaching 
their merger. Second, it is unclear what time scale should be 
used for an imminent halo merger. The dynamical friction 
time scale is relevant only aftei^the halos have started coa- 
lescing dBoylan-Kolchin et alJl2008l) . In other words, there is 
no proper definition for a "halo pair". We also note that the 
exact algorithm used for building the merger tree and the halo 
mass definition become relatively unimportant in this compar- 
ison due to the large uncertainties. Stewart et al. (2008a), who 
performed a similar analysis, have also shown that it is pos- 
sible to find a good rough agreement with observations while 
having much freedom in parameters like mass ratio and time 
scale. 

The theoretical expectation for the dark matter halo merger 
fraction brackets the observations reasonably well. This sug- 
gests that, at least within the current uncertainties, the galaxy 
merger fraction roughly follows the dark matter halo merger 
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Fig. 13. — The merger fraction as a function of redshift. Values for dark 
matter halos undergoing major mergers (with a mass ratio threshold of 3) in 
the Millennium Simulation are shown for 3 different halo masses, starling 
from the lowest halo mass we can reliably probe, that is f» 10'" ' Mq. The 
sensitivity to the chosen mass ratio threshold is indicated in the lower right 
comer A compilation of observational values that were obtained in different 
methods and with different selections nevertheless shows a rough agreement, 
also with the results from the simulation. The en'or bai's of the observations 
are largely comparable or smaller than the spread, and for the simulation 
(Poisson errors) they are very small. They have been suppressed to allow 
better readability. 



fraction. This conclusion may s eem to be different from the 
findings of i Guo & White! (|2008|) . namely that the role of ma- 
jor mergers in galaxy growth is different from that in halo 
growth, but the two conclusions are actually consistent with 
each other. iGuo & Whit3 (|2008|) found that the specific for- 
mation rate of halos through major mergers increases steeply 
with redshift but depends weakly on mass, and vice versa for 
galaxies. Nevertheless, the specific formation rate they in- 
troduced is a differ ent quantity fr om the merger fraction we 
show in Figure 

[13 i Guo & Whit3 ( 12008) defined the specific 
formation rate such that it equals the merger rate per descen- 
dant halo as defined in this work, times the Hubble time, 
with no dependence on the duration of mergers. Thus, the 
difference they found between halos and galaxies originates 
solely from the difference of the merger rate between halos 
in the dark matter simulation and the galaxies in the semi- 
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analytical model. iGuo & Whit3 (|2008|) found that the galaxy 
specific formation rate via major mergers is roughly constant 
with redshift, therefore by dividing it by the Hubble time, 
we learn that the galaxy major merger rate scales roughly as 
(1 + In comparison, they find that the halo specific for- 
mation rate via major mergers scales as (1 which after 
being devided by the Hubble time means that the halo ma- 
jor merger rate scales as (l+z)^'^- This is consistent with 
our finding that the halo major merger rate scales roughly as 
(l+zf'- However, the merger /racf/o« is proportional also to 
the merger time scale. The time scale of halo mergers T^erger 
scales roughly as (l+z)"'^ via its dependence on the Hub- 
ble time, while the observed time scale for galaxy mergers is 
probably approximate ly constant with redshift (e.g. Conselice 
120061: iKitzbichler &' White 2008). Combining these depen- 
dencies, one finds that both the halo merger fraction and the 
galaxy merger fraction scale roughly as (1 + which is in- 
deed what is seen on average in Figure [13] 

6. SUMMARY AND DISCUSSION 

We have used the Millennium Simulation to extract merger 
rates of dark matter halos. Our method differs from previous 
work in three main aspects. 

First, we reject any merger between EOF groups whose de- 
scendant subhalos, at any future time, do not belong exclu- 
sively to the same FOF group. This is done by keeping such 
FOF groups distinct until they (if at all) irrevocably merge. 
Rejecting only a fraction of such events, as in previous work, 
leads to double counting of mergers, and to false counting of 
fly-by events as mergers. Therefore, our method results in 
a lower merg e r rate , especially in the minor merger regime. 
iLudlow et al.l (l2008h find that ejections of low mass subha- 
los out of their host halos typically occur in a configuration 
where a bound group of subhalos (which was created via past 
mergers) is accreted onto a large halo, and its low mass mem- 
bers are propelled onto high energy orbits by the multiple- 
body interaction. In the context of merger counting, it is not 
clear whether the merger of the low mass halo with the group 
should be counted as a merger, or merely regarded as an inter- 
action that was interrupted by the accretion of the group onto 
the large halo. Our method assumes the latter for reasons de- 
scribed in ^2.21 Therefore, in a more conservative approach, 
our results can be regarded to as a strong lower limit for the 
merger rate. 

Second, we define each halo's mass as the mass of all the 
particles gravitationally bound to it, rather than of all the par- 
ticles constituting the FOF group. This definition reduces the 
inferred merger rate by typically 20%. The motivation for this 
definition is our finding that the total (bound & unbound) mass 
of FOF groups artificially increases by up to 50% as they ap- 
proach more massive FOF groups on their way to merge with 
them. This effect distorts the appropriate mass ratios of merg- 
ers, thereby changing t he m erger rate. A detailed discussion 
of this issue appears in ^4.31 

These two improvements make our inferred merger rate 
more co nsistent with the new predictions of Neistein & Dekell 
( l2008bl) based on the EPS model, in the sense that the func- 



tional dependencies of the merger rate on mass ratio, halo 
mass and redshift are very similar. There is a constant fac- 
tor of « 1 .5 between the EPS merger rate and the Millennium 
Simulation merger rate. In ^3.31 we provide a simple global 
fitting formula (equation (|9]l) for the merger rate per descen- 
dant halo that holds for z < 4 and all masses probed by the 
Millennium Simulation. 

Third, we also extract the merger rate per progenitor halo. 
This allows us to find a merger rate for the full range of mass 
ratios. For halos of any mass M, we can measure the rate at 
which they undergo minor mergers, major mergers and the 
rate at which they are being accreted as satellite halos onto 
more massive halos. This has significant implications for the 
redshift and mass dependencies of the major merger rate. We 
find the merger rate, in the regime where both halo masses 
are smaller than the knee of the mass function, to increase 
steeply with redshift and only slightly with mass, in a sim- 
ilar way to what was found in previous work. However, at 
high enough mass or redshift, the number of halos decreases 
exponentially and therefore the mass dependence of the major 
merger rate changes sign and starts decreasing with increasing 
mass. Also, the redshift dependence significantly weakens. In 
33.4| we provide an analytic expression for converting merger 
rates per descendant halo to merger rates per progenitor halo, 
which can be used for any theoretical merger rate. 

The two different definitions of the merger rate, per pro- 
genitor halo and per descendant halo, have different physical 
meanings. They also correspond to the two different observa- 
tional approaches towards measuring the galaxy merger rate, 
namely pair counting and morphological/kinematical identifi- 
cation. In ^we discuss the importance of this distinction and 
its relevance to observations. Finally, we find that observed 
galaxy merger fractions are consistent with the halo merger 
fraction in the Millennium Simulation within the large ob- 
servational uncertainties and uncertainties in parameters (like 
mass ratio or time scale) needed to compare the two. More 
refined comparisons still await significant improvements both 
in observational techniques and consistency and in the theo- 
retical treatment of baryonic physics. 
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APPENDIX 

A. AN ANALYTIC CONVERSION FORMULA PROM MERGER RATES PER DESCENDANT HALO TO MERGER RATES PER 

PROGENITOR HALO 

The merger rate per halo equals, in general, the number of merger events divided by the number of halos undergoing those 
mergers. Let us denote by rd{zd,M,x)dMdx the number of mergers^ that occur between Zp and Zd, whose descendant mass is 
M ± dM/2 and ratio is x ±dx/2. Similarly the number of mergers whose more massive progenitor mass is M ± dM /2 and ratio 
is X ± dx/2 are denoted as rj,{zp,M,x)dMdx, and accordingly for the less massive progenitor: rp2{Zp,M,x)dMdx. We will define 
the halo mass function is the following way: if n{z, > M) is the number of halos more massive than M at redshift z, then the 
number of halos in the interval M ± dM/2 is dM = Nh(z,M)dM. For simplicity, in the following derivation we assume the 

binary merger approximation. As shown below, the analytic formula we derive reproduces well the numerical results, therefore 
we conclude that for this matter, this approximation is not significant. 

Without loss of generality, we may assume dM 3> Mdx, therefore the progenitor masses of mergers whose descendant 
mass is M ± dM/2 and ratio is x ± dx/2 are Mj ± dMi/2 = -^(M ± dM/2) and M2 ± dM2/2 = -^(M ± dM/2). Therefore, 
rd(zd,M,x)dMdx is also exactly the number of mergers whose more massive progenitor mass is Mi ±dMi/2 and ratio isx±dx/2, 
as well as the number of mergers whose less massive progenitor mass is M2 ±dM2/2 and ratio is x± dx/2. Therefore we can 
write the equalities 

rd(Zd,M,x)dMdx = rp{zp,Mi,x)dM\dx 

rd(zd,M,x)dMdx = rp2(Zp,M2,x)dM2dx. (Al) 
The number of mergers per descendant halo, whose descendant mass is M ± dM/2 and ratio is x ± dx/2, is simply 

Rd(zd,M,x)dx = rd(zd,M,x)dMdx/Nh(zd,M)dM, (A2) 
and similarly the number of mergers per progenitor halo, whose more massive progenitor mass is Mi and ratio is x±dx/2, is 

RpiZp , Ml , x)dx = rpizp , Ml , x)dx/Nh(Zp , Mi ) 
= rd(zd,M,x)j^dx/Nh(Zp,Mi) 
= Rd(zd,M,x)Nh(Zd,M)m:dx/Ni^iZp,Mi), (A3) 

' Per unit time or redshift and possibly per unit volume. This has no significance for this derivation, as long as it is kept consistent throughout. 
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where equations dAll i and iA2\ were used in the first and second equahties, respectively. Scaling the mass from Mi to M, we 
finally arrive at 

RpiZpM.x) = RAzd, ^M,x)^Nh(zd, ^M)/N,,izp,M), (A4) 

which is the relation between the merger rate per more-massive-progenitor halo and the merger rate per descendant halo. 
In an analogous way, the merger rate per less-massive-progenitor halo is 

Rp2(Zp,M,x)=RAzd,ix+ l)M,x)(x+ l)Nh(ZdAx+ l)M)/Nh(Zp,M), (A5) 

but in order to correspond to x < 1 ratios as defined in ^2.3l we define Rp(zp,M,x < l)dx = Rp2(Zp,M, \)d^, therefore 

Rp(Zp,M,x<'^) = Rd(zd,(-+m,-)i- + l)N,,izd,i-+'^)M)x-VNh(Zp,M)- (A6) 

XXX X 

Equations (IA4l i & (IA6I 1 construct the merger rate per progenitor halo as defined in 32.31 from the merger rate per descendant 
halo, and can be applied regardless of the source of the given merger rate per descendant halo. 

In order to compare the two quantities, it is more convenient to examine equation (IA4b . because both definitions of the merger 
rate are quantified in its validity range, i.e. x > 1. The largest source of difference is the term Ni,{zd, '■^M)/Nh(Zp,M). For M 
smaller than the knee of the mass function, where it is roughly a power law with M, this term is constant with M for a given 
X, i.e. the trend of the merger rate with mass stays the same. As opposed to that, where the mass function begins to drop 
exponentially, this term decreases exponentially with M as well, and therefore so does Rp/Rd- 

For completeness we show also how the merger rate per descendant halo can be derived from the merger rate per progenitor 
halo: 

Rd{Zd,M,x)=Rp(Zp,(^)M,x)^Nh(Zp,(^)M)/N,,(zd,M) 
x+1 x+1 x+1 

M 1 x~^ M 
= Rp(Zp,—,-)—Nh(zp,(—))/Ni,izd,M), (A7) 

X+ 1 X x+1 x+ 1 

as well as how the two regimes of the merger rate per progenitor halo are related to one another: 

Ml. M 

Rp(ZpM,x)=Rp(zp,-,-)x-^Nh(Zp, -)/Nh(Zp,M)- (A8) 



